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Abstract. Bouttier, Di Francesco and Guitter introduced a method for solving certain classes of 
algebraic recurrence relations arising the context of embedded trees and map enumeration. The aim 
of this note is to apply this method to three problems. First, we discuss a general family of embedded 
binary trees, trying to unify and summarize several enumeration results for binary tree families, and 
also to add new results. Second, we discuss the family of embedded d-ary trees, embedded in the 
plane in a natural way. Third, we show that several enumeration problems concerning simple families 
of lattice paths can be solved without using the kernel method by regarding simple families of lattice 
paths as degenerated families of embedded trees. 



1. Introduction 

Several families of embedded trees have been studied in the literature. Binary trees, complete binary 
trees, several different families of planar trees and more generally simply generated tree families have 
been considered in a series of papers [8, 9, 20, 13, 3, 2, 17, 18, 12, 22, 14]: it has been showed that em- 
bedded trees naturally arise in the context of map enumeration and that properties of embedded trees 
are closely related to a random measure called Integrated Superbrownian Excursion. Combinatorial 
properties of embedded ternary trees where studied using bijections between embedded ternary trees 
and non-separable rooted planar maps [16, 10], where the authors studied a particular subclass of 
embedded ternary trees named skew ternary trees [16], or left ternary trees [10], which are embedded 
ternary trees with no node having label greater than zero. Using bijections between embedded ternary 
trees with no label greater than zero and non-separable rooted planar maps with n+1 edges they ob- 
tained amongst others an explicit result for the number of such trees of size n. Some other enumerative 
results for embedded ternary trees where derived in [19]. For the exact enumeration of embedded trees 
and related problems Bouttier, Di Francesco and Guitter [8], see also Di Francesco [13], introduced a 
new method for solving systems of recurrence relations. Bousquet- Melon [3] showed how this method 
can be used to derive deep results about the enumeration of embedded binary trees and families of 
embedded plane trees, and also about properties of the Integrated Superbrownian Excursion. The 
aim of this note is to continue the analysis of [19]. We use generating functions and the method of [8] 
to study a general family of embedded binary trees, rederiving and unifying several earlier results, 
and also the family of embedded d-ary trees. Moreover, we show that some enumeration problems 
concerning simple families of lattice paths, previously solved by Banderier and Flajolet [1] using the 
kernel method, can be treated using the method of [8]. This work is divided into three parts. The 
first part is devoted to the study of a general family of binary trees embedded in the plane, summa- 
rizing and rederiving a few of the enumerational results of [8, 3, 6]. The second part of this work is 
devoted to the study of embedded d-ary trees. The third part is devoted to the enumeration of lattice 
paths using the method of [8, 13], rederiving (and slightly extending) earlier results of Flajolet and 
Banderier [1]. Moreover, we use their method to (re-)derive other results. In particular, we derive the 
length generating function of three vicious walkers and osculating walkers, previously obtained earlier 
by Bousquet-Melou [5] using the kernel method, and Gessel. In the next section we we recall some 
properties of the family of d-ary trees and we discuss the (natural) embedding of {2d + l)-ary trees 
and (2d)-ary trees into the plane. Section 3 is devoted to a presentation of the method [8, 13] follow- 
ing the exposition of Di Francesco [13]. Throughout this work we use the notations N = {1,2,...}, 
No = {0, 1, 2, ... } and also Z = {..., -1, 0, 1, ... }. 

Key words and phrases. Embedded trees, Labeled trees. Binary Trees, Plane Trees, Height of Plane trees. Lattice 
Paths, Vicious walkers. Osculating walkers. 
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2. The family of d-ary trees 

The family of d-ary trees T, with d > 2, can be described in a recursive way, which says that a d-ary 
tree is either a leaf (an external node) or an internal node followed by d ordered ternary trees, visually 
described by the suggestive "equation" 

T ■•■ r 

T = □ + 



Here O is the symbol for an internal node and □ is the symbol for a leaf or external node. The 
generating function T{z) = X]n>o-^"^" °^ ^^'^ number of d-ary trees of size n satisfies the equation 

T{z) = 1 -I- zT\z), with r(0) = 1. (1) 

Concerning the series expansion of the generating function T(z) it is convenient consider the shifted 
series T{z) :— T{z) — 1. This corresponds to discarding external nodes (the empty tree) in the 
description above; we obtain simply generated d-ary trees T ), defined by the formal equation 

t = Ox^(t), with (^(t) (1 + t)^ (2) 

with O a node, x the cartesian product, and f{'T) the substituted structure. We refer to [21] for the 
general definition of simply generated trees. Let T„ denote the number of ternary trees of size n, and 
Tn the number of simply generated ternary trees of size n. By the formal description above (2) the 
counting series T(z) — X]n>i -^n-^" satisfies the functional equation 

f{z) = z{l + f{z)Y, f(0) = 0. (3) 
Due to the Lagrange inversion formula, see e.g. [15], the number of d-axy trees of size n is given by 
the so-called Fuss-Catalan numbers Cn4 = {d-i)n+i (tT) ' 

^"-[-"™- (d-l)n+l (J' and consequently T(.) ^ ^ J (4) 

Note that due to the definition the series T{z) and T{z) are related by T{z) = T{z) — 1. 

2.1. Embedded d-ary trees. By definition of d-ary trees each internal node with no children has 
exactly d positions to attach a new node, which are as usual called external nodes or leaves, see 
Figure 1. We embed d-ary trees in the plane by distinguishing between the cases of even and odd d, 
respectively. Equivalently, we can distinguish between (2d -I- l)-ary trees and 2d-ary trees, with d > 1. 
The root node has position zero. We recursively define the embedding of (2d + l)-ary and 2d-ary 
trees as follows. For (2d -I- l)-ary trees an internal node with label/position j E Z has exactly (2d-t- 1) 
children, being internal or external, placed at positions ±d, ±(d — 1), . . . , ±0. For 2d-ary trees an 
internal node with label/position j G Z has exactly 2d children, being internal or external, placed at 
positions ±(2d — 1), ±(2d — 3), . . . , ±1. Following [3], we call these embedding natural embedding of 
d-ary trees, because the label a node is its abscissa in the natural integer embedding of the tree. 




-2-1012 -2-1012 -3-2-10123 -2-1012 



Figure 1: Size one naturally embedded binary, ternary, quaternary and quinary trees together with 
their external nodes, i.e. the possible increments i E {±1}, i E {0, ±1}, i E {±1,±3} and i E 
{0,±1,±2}. 



In this note we are interested in the number of embedded d-ary trees having no label greater than 
j, with j E N. Let Tj^2d+i{z) and Tj^2d{z) denote the generating function of embedded (2d -I- l)-ary 
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and 2d-ary trees having no label greater than j, j G N, with initial values T_i 2d+i = T^2,2d+i = 
■■ ■ = T-dM+i = 1 and T_i_2d+i = T^2M+i = • ■ ■ = T_2d+i,2d+i = 1- Following the observation of 
Bousquet-Melou we can think of Tj(z) as the generating function of embedded d-ary trees with root 
labeled j. By definition we obtain the following system of recurrences for Tj{z)^. For {2d+ l)-ary trees 
we get the system of recurrences 

d 

T,{z) - 1 + 2 n ^^(^)' ^ ^ ^' ^^th T_,(z) = 1, for l<j<d, (5) 
e=-d 

and for 2d-ary trees we get the system of recurrences 

d 

T,{z)^l + zl[(T2e-i{z)T^2e+i{z)), j>0, with r_,(z) = 1, forl<i<2d-l. (6) 

£=1 

Note that for both cases we have 

Tj{z) T{z) for j oo, 

in the sense of formal power series, where T denotes the overall generating function (1) of {2d+ l)-ary 
and 2d-ary trees, respectively. Note that this observation turns out to be crucial for the solution of the 
recurrence relation; see the original paper of Bouttier et al. [8] and the next section. In the work [19] 
a different embedding for 2(i-ary trees is suggested. However, the embedding above for 2d-ary trees 
turns out to be more easily analyzed and more natural, since the nodes are evenly placed in the plane. 

3. A METHOD FOR SOLVING INFINITE SYSTEMS OF ALGEBRAIC RECURRENCE RELATION 

Bouttier, Di Francesco and Guitter introduced a method for solving certain classes of algebraic re- 
currence relations arising the context of embedded trees and map enumeration. Our presentation of 
their method follows the exposition of Di Francesco [13]. For a given integer A: G Z let Tj{z), with 
j > k, denote a family of generating functions. Assume that the Tj{z) satisfy algebraic recurrence 
relations expressing Tj{z) in terms of a finite number of previous terms Tj-i{z),Tj-2{z), . ■ . , Tj^d{z), 
with d G N. The boundary data needed to entirely determine Tj{z) should consist of d consecutive 
initial values of Tj{z). Assume further that in the sense of formal power series limj^oo Tj(z) exists, 
with limj^oo = T{z); note that T(z) is also the solution of the unrestricted recurrence rela- 

tion for Tj{z), holding for all j G Z. Exploiting the fact that liuij^oo Tj{z) = T{z) one uses the 
ansatz Tj{z) = T(z)(l — pj(z)), where Pj{z) denotes an a priori unknown formal power series with 
liiUj^ao Pj{z) = 0. This allows to linearize the recurrence relations at large j, similar to first order 
asymptotic series expansion. 

A first order expansion of the recurrence relation for Tj(z) in terms of Pj{z) leads to linear recurrence 
relations for Pj{z) — p^p{z). It is readily solved using the classical ansatz pj^'(z) = a ■ , with 
unspecified a. We can deduce that the general solution of the linearized recurrence relation is given 
by pf\z) = Yl'e.^i '^i ' -^£1 where the Xi{z), with 1 < £ < d, are all solutions with modulus less one of 
the characteristic equation of the linear recurrence relation for the first order approximation pj (z) = 
p'^p [z). In order to obtain the solution of the original problem one uses a full asymptotic series expan- 
sion of the recurrence relation for Tj{z) in terms of pj{z) = '}2m,...,nj^>Q "ni,...,nd"ni, Ilfci {^lY'^ 
and compares order by order the contributions to the true solution. We recursively obtain the unspec- 
ified coefficients ani....,nfc, usually depending on Xi(z), 1 < i < d, with free parameters ae^, where ei 
denotes the ^-th unit vector. 

The main difficulty is solve the recurrence relation for the coefficients — ctm,. • Once these re- 
currence relations are solved, one can hopefully derive a compact expression for Pj{z) and subsequently 
adapt the unspecified parameters ae^, 1 < ^ < d, to the initial conditions Tj-i(z), Tj-2{z), . . . , Tj-d{z). 



^Subsequently, we will usually drop the subscripts 2d+ 1 and 2d of Tj^2d+i{z) and Tj 2d(^) iu order to simplify the 
presentation. 
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4. General families of embedded binary trees 

The family of ordinary (incomplete) binary trees 7i, enumerated by the Catalan numbers, whose 
counting series T = T{z) ~ X^tsTi satisfies the functional equation 

T{z)^l + zT{zf. 

Bousquet- Melon [3] considered the embedding of this tree family in the plane according to 

Tj{z) = 1 + zrj_i(z)r,+i(z), J e z. 

Here Tj{z) denotes the generating function of a tree with root at position j £ Z. Bouttier et al. [8, 9] 
and Bousquet-Melou [3] studied two families T2 and of embedded plane trees which are closely 
related to families of maps. They can be realised as certain families of embedded binary trees. Let 
F = F{z) = X^TsTg ^^"^ ^ ~ G{z) — J^TeTs denote the counting series of the families T2 and 
7^, satisfying the functional equations 



l-2zF{z)' ' ' l-3zG{z)' 
or equivalently, 

F{z) ^l + 2zF{zf, G{z) = l + 3zG{zf. 
These tree families are embedded according to 

^'^'^^ l-z{F,.,{z)+F,^,{z))' ^'^'^ ^ 1 - z(G,_i(z) + G{z) + G,+i(z)) ^ ^' 
or equivalently by 

F,{z) = 1 + zF,{z){F,^^{z) + G,(z) = 1 + zG,(z)(G,_i(z) + G,(z) + G,+i(z)). 

For the three tree families 7i, T2 and T3 it was shown that the generating functions of trees with 
small labels, i.e. tree in which all labels are less or equal j, are algebraic and explicit expressions were 
obtained. 



4.1. Embedding of a general family of binary trees. We discuss properties of the family T 
of weighted binary trees, defined according to a functional equation for its counting series T — 

T{z,Vi,V2,Wi,W2,W3) = EggT ^'"^'i 

T = 1 + z{2vi + V2)T + z{wi + W2 + 2w3)T^. 

We can interpret f2, wi, 1/72, W3 either as weights, vi,V2,wi, 102,1113 > 0, or as variables encoding 
different kinds of nodes, which would lead to a refined enumeration of trees. Concerning the second 
point of view one could for example consider [z^v^^v^^Wi^w^w^^]T, with mi + m2 + ^1+^2 + ^3 — n. 
By solving the quadratic equation for T one easily obtains the following explicit result. 

1 - z{2vi + V2) - ^(1 - z{2vi + V2)Y ~ 4:Z{wi +W2+ 2W3) 
2z{wi +W2+ 2W3) 

We reobtain the previously considered families and several other tree families, binary and non-binary, 
by suitable sometimes non-unique choices of vi, V2 and Wi, W2, W3. 

Example 1. Binary trees (Catalan numbers) A000108 are obtained by setting vi = V2 = W2 = 
W3 = and wi — I, the number of rooted Eulerian edge maps in the plane A052701 are obtained 
by setting vi = V2 — wi = W2 — and W3 = 1, Blossom trees or equivalently rooted planar maps 
A005159 are obtained by setting vi — V2 — wi — and ti;2 = W3 = 1, Schroder trees (large Schroder 
numbers) A006318 can be obtained setting V2 — W2 — W3 — and vi = wi — I, planar rooted trees 
with tricolored end nodes A047891 can be obtained setting vi — W2 = W3 — and V2 — wi = 1, the 
choice vi ^ V2 = wi ^ 1 and W2 — W3 — gives sequence A082298, the choice vi = W3 = 1 and 
V2 = wi — W2 — gives sequence A103210; several other sequences in Sloane's Encyclopedia [23] can 
be obtained by suitable choices of the parameters. 
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We embed this family according to the following recurrence relation for Tj ~ Tj{z,vi,V2,Wi,W2,W3). 

Tj = 1 + z(uiTj_i + viTj+i + V2Tj) + z(^wiTj^iTj+i + W2Tf + wsTj (Tj_i + T^+i)^ , (8) 

with j e Z. We will see that we can reobtain the previously discussed families 7^, 7^ and and their 
counting series by the following choices of the weights/ variables Wi, W2, Wi,W2, w^: T{z, 0, 0, Wi, 0, 0), 
T{z, 0, 0, 0, W2, W2) and T(z, 0, 0, 0, 0, ^2)- 

We will show that for several choices of the weights Wj and arbitrary weights Vi the generating functions 
of trees with small labels in the embedded family T, i.e. tree in which all labels are less or equal j, 
can be explicitly obtained. 

4.2. Trees with small labels. Our starting point is the recurrence relation below for Tj. 

= 1 + z(uiTj_i + viTj+i + V2Tj) + z{wiTj^iTj+i + W2Tf + w^Tj (Tj_i + T^+i)^ , (9) 

for j > with initial value given by T_i = 1 or r_i = 0, depending on particular counting problem, 
see [8, 9, 13, 3]. Following the approach presented in Section 3 we use that fact that for j tending to 
infinity we have Tj T in the sense of formal power series, with T given by (7). We make the ansatz 
Tj = r(l — pj), where T = T(z, ui, W2, wi, 102, W3) denotes the generating function of the family T 
defined by (7), with pj ^ as j tends to infinity. We expend Equation 9 with respect to the ansatz 
and compare the terms tending at a similar rate to zero in the asymptotic expansion of Tj as j tends 
to infinity, neglecting terms pjPj^i and pjPj-i. We get the linearized equation 

-Tpj = -zT{vi{pj-i + pj+i) + V2Pj) - zr^(wi(pj_i + pj+i) + 2w2Pj + wz{pj-i + 2 + pj + p^+i)). 

Now we make a refined ansatz pj — in order to solve this linear recurrence relation for pj , assuming 
that X is a formal power series depending on variables/ weights z,vi,V2,wi,'W2,W3 with \X\ < 1. We 
obtain the so-called characteristic equation for the series X, 

l = z(yi{^+X)+V2)+zT(wii^+X)+2w2 + W3{^ + 2 + X)y (10) 

We observe that X is a power series in z, wi, U2, wi, 102, W3 and has non-negative coefficients. Conse- 
quently, the proper solution is given by 

1 - Z{V2 + 2T{W2 + W3)) - y^(l - Z{V2 + 2T{W2 + Wi))f - 4z2(i;i + T{wi + W3))2 

X = -— . (n) 

2z[vi + T[wi + W3)) 

One readily checks that the expression above for X is indeed a power series in z, vi,V2,wi,W2, and 
has non-negative coefficients. Using the definition of the series T we can express T solely in terms of 
the series X 

ti{X) + ^ti{XY + U2{X){vi{l + X^) + V2X) 

h{x) ' 

with respect to the polynomials ti{X) — ti(vi,wi,W2,W3, X) and t2{X) — 12(^01, W2, w^, X) defined 

by 

ti{X) = wi{l+X^)+2w2X+W3{l+Xf^vi{l~Xf, t2{X) = wi{l~X+X^)+W2X+W3{l+X^). 

We make the more refined ansatz pj — X]i>i with unspecified ai and — ai{X), which 
amounts to an asymptotic expansion of pj for j tending to infinity. Next we compare the terms with 
the same order of magnitude in (9) as j tends infinity. We obtain from (9), using the relation (10), 
the following recurrence relation for an+i, with n > 0. 

i=l 
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We observe that the variable V2 only appears in the defining equations for series T and X, but not 
in the recurrence relation for q;„. Introducing the quantity /3„+i = q;„+i(^ + wi + w^j , n > 0, we 
obtain the simplified recurrence relation 

1 1 " 1 

Pn+i + X^+' - ^ - ^) - E + ^2 + (— + X')) . (13) 

i=l 

Let f{t) denote the formal power series f{t) — X]n>i Pnt^- Equation 13 is equivalent to a functional 
equation for f(t): 

f{tx) + /( A) - (1 + x)f{t) = w,f{tx) •/(!)+ w2f{tY + w,f{t) ■ [f{tx) + /(I) 

One already knows the solutions of Equation 13 in the cases z«2 = W3 = 0, see Bousquet-Melou [3], and 
Wi = 0, Wi = W2 — 0, see Bouttier et al. [8] and also [3]. We will provide the solution of Equations 12 
and 12, respectively, in the case wi and W2 — W3, excluding the degenerate case wi — W2 — wj, — 0. 

Lemma 1. For given parameters wi and W2 = w^, excluding the degenerate case Wi = W2 — = 0, 
the solution Q!„ of the recurrence relation 12 is for n> 1 given by 

_ X''-^a'i{wiX + W2{1 + X + X^))"''(l - X") 

+ + W2)"'\l - X)2«-l(l +X + X2)"-l(l + ' 

We could not solve directly the functional equation for f(t). Instead we obtained the solution in an 
experimental way using the computer algebra software Maple. Once the solution of the recurrence 
relation is guessed, it is readily rigourously checked that it satisfies the recurrence relation (12), or 
equivalently that the generating function f{t) — J2n>i f^^t" satisfies the stated functional equation. 
Unfortunately, we could not solve the recurrence relation in full generality W2 7^ w^, except for the 
already known special case wi = u'2 = and 7^ [8, 3]; it is given by 



'(1 _ X)2"-i(l + X + X2)«-i(l + X) 



However, the result of Lemma 1 already covers and generalizes the result for two previously treated 
families, the cases ui = W2 = W2 = of binary trees and ui = f2 = wi = of a family of planar 
trees, which we interpret as embedded binary trees. It seems that the structure of the values a„ is 
not regular in the other cases. We performed some computer experiments and we state the following 
conjecture on the values of a„ for ^2 = and wi — W3 = 1. 

Conjecture 1. In the case W2 = and wi — — 1 the solution an of the recurrence relation 12 is 
given by 

a^X"-^Pn{X) 

+ 2)""^(1-X)2«-2(1+X)4^J(1+X2)L'^J ' 

where the sequence of polynomials (p„(X))„gN with initial values 

Pi{X)^l, P2{X) = l, P3(X) =X4 + 2X3 + 2X + 1, 
is for n > 2 recursively defined by 

P2n{X) = P2n-l{X) - 2X^p2n-2{X), P2n+l{X) - p„+2 (X)p„+i (X) - 4X V (^)P«-1 (^) ■ 

We return to our previous case of wi and W2 — wj,. In order to simplify the presentation we set 

(^+^l+^2)(l-X2)(l-X3) 

"1 = 7 ^7-^ 1^ , V , V2\\ \ = X{X,Vi,Wi,W2), 

(wiX + W2{1 + X + A2) ) 

and obtain the following result. 
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Theorem 1. For given parameters Wi and W2 — w^, excluding the degenerate case Wi = W2 — — 0, 
a solution of the recurrence relation 8 with free parameter A is given by 



Tj = T ■ [1- 



{wiX + W2(l +X + X2))(l - AXJ + i)(l - / ' 



with series X given by (11). 



Now we can easily reobtain the previous results of [8, 3] by suitable choices of Vi,V2,Wi,'W2-,W3 and 
adapting A to the initial value T^i- The quadratic equation relating A and T_i normally has two 
distinct solutions; we use the fact that Tj a priori has a power series expansion at 2; = to identify 
the right solution. 

Corollary 1 ([8, 3]). In the case of embedded binary trees, Vi = V2 = W2 = = and Wi = I with 
T_i — 1, we reobtain the result 

Tj = T- ): z../.): i>-i. 



(1 


-XJ+2)(1 


- 


(1 


-XJ+4)(1 






= V2 = Wi 


= and 


(1 


-XJ+i)(l 




(1 


-XJ+2)(1 





the result 

Tj^T-^ 'i-^Si ^77T^' j>-l- 

Remark 1. As mentioned above one can readily obtain numerous enumerative results from The- 
orem 1. The solutions turn out to be usually more involved due to the adaption to initial values 
T_i = 1 or r_i = 0. 

4.3. The height of planar trees. A more general form of recurrence relation (9) reads the following 
way. 



Tj = 1 + z(wiTj_i + V2Tj + vsTj+i) + zywiTj^iTj+i + W2Tf + w^TjTj^i + WiTjTj+i] 

It seems very difficult to obtain solutions for this recurrence relation. However, there exist a subclass 
V3 = wi = W2 — Wi — 0, setting for the sake of simplicity W3 — 1, which is explicitly solvable 

Tj = l + z{viTj^i+V2Tj)+zTj^iTj, for j > 1. (14) 

This was observed earlier by Bousquet-Melou [6]. This subclass is of particular importance due to 
the connection with the height of plane trees [11], corresponding to the case vi — V2 — 0, with initial 
condition Tq{z) — 1. By the approach presented in Section 3 we use the fact that 

Tj^T, with T ^1 + z{vi+V2)T + zT^, 

and the ansatz Tj = T(l — pj). This leads to 

-pjT = -zT(wi/9j_i + V2Pj) - zT^(pj^i + pj) + zT^pj-ipj. 

For the first order expansion we consider the terms tending at the same rate to zero as j tends to 
infinity, neglecting the term pj^ipj. We get the linearized equation 

Pj = z{vipj^i + V2pj) + zT(pj^i + Pj). 

This recurrence relation is readily solved by pj = X^ , with X = X{z) given by satisfying 

z{vi + T) 



X = 



1 - z{v2 + T) ■ 



As before he more refined ansatz pj — X]i>i cti{X^y , with unspecified ai and = ai{X) leads to a 
recurrence relation for 



X" ^ 1 

, n>0. 
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Proceeding as before, we introduce the quantity /3„ = a„(^ + n > 1 and solve the arising 

recurrence relation in an experimental way using the computer algebra software Maple. We obtain 
the solution 



n > 1, 



with unspecified ai. In order to simplify the presentation we set ai = A(l — + 1) and obtain 

the following result. 

Theorem 2 ([6]). A solution of the recurrence relation 14 with free parameter X is given by 

/ (f + 1)(1-X)AX^\ 
' \ (l-AX^+i) > 

with X = z{vi + T)/(l - z{v2 + T)). 

Setting vi ^ V2 = and adapting to the initial condition Tq = I gives the following result. 
Corollary 2 ([11, 6]). The generating function of plane trees of height <j is given by 

1 - X-J'+i 

r,=r-^-— ^, with X = zT/{l-zT)=T-l. 
Note that the free parameter A in Theorem 2 allows to obtain refined enumerations, as shown in [6] . 

4.4. A family of ternary trees. We have seen that one can eliminate the variables vi and V2 
from the recurrence relation (12) (and also in (14)) for by a proper substitution, leading to the 
simplified recurrence relation (13) for the values /?„. This is not the case anymore even for ternary 
trees. Consider for example the family T of weighted ternary trees, defined according to a functional 
equation for its counting series T = T{z, vi,V2) = X^tsT'^''^'' 

T{z) = 1 + z{2vi + V2)T{z) + zT{z)\ 

embedded according to recurrence relation 

Tjiz) = 1 + z{viTj^i + viTj+i{z) + V2Tj{z)) + zTj^i{z)Tj{z)Tj+i{z), 

with j G 1. Proceeding as before, i.e. making the ansatz Tj — T(l — pj) and subsequent refinements 
Pj = J2n>i (^'ri{X^y^ with X being the solution of 



l^z{v,{^+X)+V2)+zT'(^^ 



zT'{- + l + X 



with \X\ < 1, one obtains the recurrence relation 



1 1 " 1 



il+j2+i3=" 



Unfortunately, we are not able to solve this recurrence relation for vi ^ 0. We observe that as before 
the variable V2 only appears in the defining equations for series T and X, but not in the recurrence 
relation for a„. In the case wi = we can use the solution of [19], and subsequently may obtain a 
refinement of a result of [19]. 
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5. Embedded (2d+1)-ary trees with small labels 

The starting point of our considerations is the system of recurrences (5). We make the ansatz Tj(z) = 
T(l — pj), with pj = Pj{z) — > as J tends to infinity, for z near zero. We expend Equation 5 with 
respect to the ansatz and obtain 



By definition of (2(i+ l)-ary trees (1) we have 1 = T — zT^'^^^. Consequently, 

d 

r(i - p,) = r(i - zT^" + zT^' n (1 - P^^^)) ■ 

i=-d 

The equation above is equivalent to 

d 

l-p^^l-zT^'i + zT^'^ J] {l-p,+,). 

l=-d 

By expansion of the product on the right hand side of the equation above we obtain the main equation 

d 2d+l , I s 

p, = zT^' - zT^' n (1 - ^.+^) = E (-1)'" E ( n p^^b^ ) ' (15) 

with — {bi, . . . , be} running over all subset of {— rf, . . . , d} of size £, 1 < £ < 2d+ I. Comparing the 
terms tending at a similar rate to zero as j tends to infinity we obtain the linear recurrence relation 

d 

Pj = ^^^"^ E P:>+^- 
i=-d 

An ansatz pj = aX^ , assuming that there exists a formal power series X — X{z) with \X\ < 1 for z 
near 0, leads to the so-called characteristic equation 

d d 

1 = zT^'^ X^ or equivalently 1 = Z ^ with Z := zT^"^. (16) 

e=-d e=-d 

The equation above is identical to the characteristic equation of lattice path with step set S = {{!,£) \ 
—d < £ < d}, see Banderier and Flajolet [1]. We can use the very general considerations of [1] 
summarized below in Lemma 2 concerning such equations. 

Lemma 2 (Banderier and Flajolet [1]). Let 5 C Z denote a non-empty finite subset of the integers 
S — and H := {wi,...,Wm} Q the set of associated weights. The characteristic 

polynomial P{X) associated to S and 11 zs a Laurent polynomial in X given by P{X) — X^fci w^X^'^ . 
Let c = — minj^f} and d = maxj^^}. The characteristic equation associated to S is given by 

1 - ZP{X) = 0, or equivalently X' - ZX^PiX) = 0. 

For Z near zero the characteristic equation has c + d solutions, of which c solutions Xi(Z), . . . Xc{Z) 
are small solution with \Xf(Z)\ < 1 for Z near zero. These c so-called small branches are conjugate 
of each other at Z = 0: there exist two functions A and B analytic at Z = and nonzero there such 
that in a neighbourhood of zero one has 

Xe = Xi{Z)=Lo^-^Z^/''A{Lu'^~^Z^^''), with = e^"/^ 

where i denote the imaginary unit — —1. 

The result of Banderier and Flajolet [1] was initially derived in the context of the enumeration of 
(weighted) lattice paths. We apply Lemma 2 to case S — {£ \ ~d < £ < d], wg = 1, —d < £ < d 
and Z = zT'^'^. Consequently the equation 1 = ZY^'l^_^X^ , with Z = zT"^"^, has d small solution 
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Xi{z), . . . Xii{z). We refine the previous ansatz pj — aX^ in terms of the d solutions Xi{z), . . . Xii{z) 
of the characteristic equation in the following way. 

ni,...,nci>0 n>0 

with ao,o,...,o = and unspecified initial values agj,, 1 < ^ < d, where eg denotes the £-th unit vector. 
By the refined ansatz the main equation (15) reads the following way. 

2d+l I , V 

E «nX^" - 7^T^^ E (-1)'"' E n ( E «nX(^+'"=)" j . 

n>0 l=\ bfC{-d,...,d} A;=l ^n>0 ' 

Our goal is to determine the unknown coefficients a-a — an(X) as functions of X\, . . . ,Xd and the 
unspecified initial values ae,, , 1 < i < d. In order to do so we compare the terms with the same order 
of magnitude in (9) as j tends infinity; this corresponds to some kind of coefficient extraction with 
respect to [X^*^"^]. We obtain the for n = (ni, . . . ,nd) € Nq \ {0,ei, . . . ,6^} the recurrence relation 



d 2d+l 

/ 1 

a 



here e Nq \ {0} denotes a vector of length d, for 1 < fc < ^ and 2 < i < 2d + 1, hi ^ {bi, . . . ,b(} 
runs over all subset of {—d, . . . ,d} of size £, 2 < £ < 2d + 1. Let /(w) denote the formal power series 
/(w) = X]n>o "^n^"- recurrence relation above for a„ = an(X) is equivalent to a functional 

equation for /(w) 



2d+l 



L/(w) + E /(X'^) = E (-1)' E ( n /(X'^w) ) ; (19) 



^J^2d • 

i=-d e=2 b,C{-d,...,d} ^fe=l 

here we use the notation X^w :— {Xfwi, X2W2, ■ ■ ■ , X^Wd)- A crucial step towards solving recurrence 
relation (18) is to study the recurrence relation a„ with n = rifee^. = (0, . . . , 0, rife, 0, . . . , 0) for 1 < 
k < d. Note that the simplified recurrence relation involves only terms Q!„ such that ni = for £ ^ k. 
Once these one parameter solutions are obtained, the general solution is immediately determined 
by equation (18). Moreover, by definition the formal power series /3j(X) is a solution of the main 
equation (15). 

5.1. One parameter solution. A solution with only one free parameter can be obtained by using a 
simpler ansatz using only one of the series Xi{z), . . . ,Xd{z). Equivalently, we consider an = an(X), 
with n = (0, . . . , 0, rife, 0, . . . , 0) for 1 < k < d. We obtain the following simple solution. 

Lemma 3. Let X = X{z) denote any of the d small solutions Xi(z), . . . , Xd{z) of the characteristic 
equation 1 = -zT^'' X^. The recurrence relation (5) admits a solution with free parameter A 

^ (l-AX-'+i+^)(l-AX^'^+3+^) 
^ {I - XXd+2+i){l - XX^d+2+jy • 

Proof. Since by definition the series X satisfies l/izT^'^) = t^^e recurrence relation (18) 

simplifies to 

d d 2d+l / ^ \ / ^ \ 

a«+i(-E^^+E^^^""^0 = E(-i)^ E (n«.o ^ (n^^^^O' -^1- 

e=-d e=-d 1=2 Efc=i9fc="+i '^=1 b«c{-<i,...,rf} ^fc=i ^ 

By experiments with Maple we obtain a solution with free parameter ai, 

<X"-i(l-X"'') 
^ (1 - X^)(\ - X)"-i(l " X'i+i)"-! ' ^ ^ 
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One can checks for small d = 1,2,3,... that the arising function f{w) — J2n>i'^nw"' satisfies the 
functional equation 

d d 2d+l , t 

-f{w) E + E /(^'-) - E (-1)' E ( n /(^'^-) 

i=~d l=-d 1=2 b4C{-d,...,d} ^fc=l 

Now we set a\ = AX'*+^(1 — X)(\ — X'^^^') in order to simplify the calculations. We get T,- = 
T(l — X]n>i '^"^"'"i o'^"^ readily checks that the stated solution satisfies recurrence relation (5). □ 

5.2. The general solution. An immediate application of Lemma 3 and the explicit result (20) for 
a„ is the following result. 

Proposition 1. Let pj — X]n>o Q^nX-'", with pj = pj{z) ~ pj(X.), and coefficients — Q;n(X) given 
by 



1 



/orn = (0,0, ...,0), 



1<'^r" (1 - ) 
; ; for n = ntOi, for ng > 1, 1 < £ < d, 
determined by recurrence relation (18) for n G Nq \ {0,ei, . . . ,ed}, 

with unspecified initial values aee, I < i < d. Then, the formal power series Tj — T{\ — pj) satisfies 
the recurrence relation 

d 

T,{z) = l + z n T,+,{z). 

e=-d 

Equivalently, the formal power series pj satisfies the equation 

2d+l , t 

e=l beC:{-d,...M} ^ k=l 

Proof. By the recursive description (18) of and the simple observation that for z near zero we have 

1 k 

= -zT^dx'''^ y (zT^'^ y x'A , 

zT^d ~f l^i=^d^ k>0 £=0 

the values an = can be written as formal power series in T = T{z) = T — 1, according to 
zT'^d. = X - T = - f, with f{0) = 0, and X. Consequently, by definition of the values (18) the left 
and right hand side (15) coincide. □ 

The huge obstacle concerning our enumeration problem (5) is to explicitly determine the values in 
the general case in order to adapt the initial values A^, 1 < I' < d, to the initial conditions T_£ = 1, 
1 < £ < d. The only way known to us to obtain an is either guessing the solution after experiments, 
or to solve the functional equation (19). Unfortunately, we do not know how to directly solve (19) 
and we did not manage yet to guess a general formula for an- 



6. Embedded 2d-ary trees with small labels 

The considerations for embedded 2(i-ary trees are similar to (2d + l)-ary trees, therefore we will be 
more brief. According to the ansatz Tj — T{1 — pj) we expend Equation 5 and obtain the equation 

d 2d , I ^ 

= ^zT^<^-^ + zT^'^-i - P.+2.-i)(l - P,-2£+i) = zT^'^-^ Y^{~lY ^ n Pi+^>' ) ' (21) 

i=l e=l bi<ZB2d ^k = l ' 
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with hi ~ {bi, . . . , bg} running over all subset of — {~(2d — 1), — (2(i ~ 3), . . . ,2d — 1} of size £, 
1 < i < 2d. Comparing the terms tending at a similar rate to zero as j tends to infinity we obtain 
the linear recurrence relation 



An ansatz pj = aX^ , assuming that there exists a formal power series X = X{z) with \X\ < 1 for z 
near 0, leads to the characteristic equation 

d d 

1 = zT^'^-i {X^''^+X-^'+^), or equivalently 1 = {X^'-^ + X-^'+^) , with Z := zT^''-\ 

1=1 i=i 

(22) 

We apply Lemma 2 to case S = {2£ ~ 1, — 2£ + 1 \ 1 < £ < d}, with weights all equal to one, and 
Z = zT'^'^-^. Consequently the equation 1 = ZY.j^^ {X^^^'^ + X-'^'^+^) , with Z = zT^-^-i, has 2d- 1 
small solution Xi{z), . . . X2d-i{z). As before, we refine the previous ansatz pj — aX^ in terms of the 
2c? — 1 solutions Xi(z), . . . X2d-i{z) of the characteristic equation in the following way. 

Pj — 2^ "ni,...,n2d~l^l ■■■^2d-l — / ^ , 

ni,...,n2d-i>0 n>0 

with ao.o,...,o = Sind unspecified initial values Ofe^, 1 < ^ < 2d — 1, where denotes the £-th unit 
vector. According to the refined ansatz the equation (21) we obtain the for n = (ni, . . . ,n2d-i) G 
Ng''"^ \ {OjGi, . . . ,e2d-i} the recurrence relation 



a, 



^ a 2d / ^ \ / *^ 

— 1 D — O -sr^P V 7 1 / U f?_ . V J 1 



(23) 

here e No'^"^\{0} denotes a vector of length 2d-l, for 1 < fc < £ and 2 < ^ < 2d, = {5i,. . . ,bt\ 
runs over all subset of Bad = {-(2d - 1), -(2d - 3), . . . , 2d - 1} of size 1 < £ < 2d. 

6.1. One parameter solution. The solutions with one free parameter can be obtained by using a 
simpler ansatz using only one of the series Xi(z), . . . , X2d-\{z)- Equivalently, we consider «„ = an(X), 
with n = (0, . . . , 0, rifc, 0, . . . , 0) for 1 < A: < 2d — 1. We obtain the following simple solution. 

Lemma 4. Let X = X{z) denote any of the 2d — I small solutions Xi{z), . . . , X2d~i{z) of the 
characteristic equation 1 — zT^'''^^ '^'j^^{X^^^^ + X^^^^^) . The recurrence relation (6) admits a 
solution with free parameter X 

(1 - XX'^+^+^){l ~ AX^'^+^+J) 

■ (1 - XXd+3+j)(l - AX3'i+2+J) ' • 

Proof. By experiments with Maple we obtain a solution with free parameter ai of the simplified 
recurrence relation (23) 



j,nj^2(n-l)(^-l^ _ jj^n(2d-l)-) 
(1 - X2'i-l)(l-X2)"-l(l-X2'i+l) 



rt > 1 (2A) 

2d-lUl _ Y2\n-l(1 _ Y2d+l\n~l' - ' ' 



Now we set ai = XX'^^^{1 — X'^){1 — x'^'^^^) in order to simplify the calculations; one readily checks 
that the stated solution satisfies recurrence relation (6). □ 

6.2. The general solution. An immediate application of Lemma 4 and the explicit result (24) for 
ttn is the following result. 
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Proposition 2. Let pj — J2n>o c^nX^"; with pj 
by 



Oil, 







2{ne-l) 



nt{2d-l) 



(1 



Pj{z) — pj(X), and coefficients an = ctnO^) given 

/orn= (0,0,...,0), 

— for n = n^e£, /or ng > 1, 1 < £ < 2d — 1, 



determined by recurrence relation (23) 

with unspecified initial values ae^. , 1 < i < 2d - 
satisfies the recurrence relation 



/orneN2'^-i\{0,ei,...,e2d-i}, 
1. Then, the formal power series Tj — T(l — pj) 



7. Enumeration of lattice paths and degenerated embedded trees 



Fix a set of step vectors S — {(1, bi), . . . , (1, 6m)} 
called a walk, is a sequence (wi, . . . , «„) such that 
restricted to No x Nq. An excursion is a meander 
is often useful to consider weighted lattice paths 
to step (1, 6f). The weight of a path is defined as 



with 6£GZ, l<^<r7i. A simple lattice path, also 
for each Vi € S. A meander is a simple lattice path 
with starting point and end point on the y-axis. It 
n — {wi, . . . , Wm}, where weights wi is associated 
the product of the weight of the steps. 



pa t h/ wa I k 



(constrai ned) 















_LJ LJ LJ L 












. 1 L /j _ ,0 1 ^ ^ / _ - li^ - 


"1 r n pL -\ r 
























me a n d e r 


e X c u r s i on 



constrai ned) 



Figure 2: Degenerated 1-ary tree defined by T{z) = 1 + zT{z), and three types of lattice paths: 
unrestricted paths, meanders restricted to Nq x Nq, and excursions starting and ending at level zero. 



Banderier and Flajolet [1], amongst many other things, derived the generating functions of meanders 
and excursions with respect to a set of step vectors S and weights H using the kernel method. In 
the following we will rederive (and slightly refine) these generating functions using the method of 
Section 3. Following [1] we introduce the characteristic polynomial P{X) := X)fci WiX^^ of step set 
S = {(1, 6i), . . . , (1, 6m)}, with c = — min{6£} and d = max{6£}, and weights 11 = {wi, . . . , Wm}- Let 
Tj{z) denote the generating functions of all weighted meanders with steps S = {(1, 6i), . . . , (1, 6m)} 
and weights 11 — {wi , . . . , Wm} starting at level j, with j > 0. We have the infinite system of recurrence 
relations 

m 

Tj(z) = 1 + zY,weTj+be{^), 3 > 0, (25) 

with initial conditions r_i(z) = • . . = T^c{z) = 0. This recurrence relation can be interpreted as a 
recurrence relation for embedded trees with respect to the class of degenerated unary trees defined by 
the equation T(z) = 1 + zT{z). 

For j — !■ oo we have convergence in the sense of formal power series Tj{z) T{z), where T(z) is the 
generating function of unconstrained lattice paths, i.e. starting at zero and ending anywhere, with 
steps S = {(1, 6i), . . . , (1, 6m)} and weights 11 = {wi, . . . ,Wm}, 

m ^ 

r(z) = l + z^w,r(z), or Tiz)= _ 

1=1 ^ ^^^'^^ 
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where P{X) = J^JLi w^X^^ denotes the characteristic polynomial of steps and weights. We use the 
ansatz Tj{z) = T(z)(l — Pj{z)) to obtain 



- p,{z)) - 1 + zY,^,T{z){l - 
1=1 

Consequently, we get a linear recurrence relation for Pj{z), 

m 

Setting Pj{z) = X^ we obtain after simple manipulations the characteristic equation 

m 

1 - z ^ wiX''' = 0, or equivalently, 1 - zP{X) = 0. (26) 



1=1 



By Lemma 2 there exist c solution small solutions Xi(z), . . . , Xc{z) for z in a neighborhood of zero, 
and the general solution is given by 

c 

Pji^) ^^OiiXl, 

1=1 

with unspecified a^, 1 < ^ < c. Consequently, we obtain the following result. 

Theorem 3. A family of solutions of the system of recurrence relations (25), well defined for z in 
a neighbourhood of zero, is given in terms of the c small solutions Xi, . . . ,Xc of the characteristic 
equation (26), 

c 

T,{z)^T{z){l-Y,o^,Xl), 

1=1 

with free parameters ai, with I < £ < c. The parameters ai are independent of j , but may depend on 
Xi, . . . , Xc and z. 

Adapting to the initial conditions of the meanders r_i(z) — ■ ■ ■ = T_c{z) — lead to a system of c 
linear equations 

c 

1 - ^ aiX][^ = 0, for l<i<c. 

1=1 

This system is easily solved using Cramer's rule and Vandermonde's determinant; we obtain the result 

]\l<i<k<ci^i^ ^k)) 

ni<i<fe<c(^* ^ 

Consequently, we get after simple manipulations the following result. 

Corollary 3. The generating function of meanders with steps S = {(1, 6i), . . . , (1, 6r„)}, weights 
n — {t«i, . . . , Wm} and characteristic polynomial P{X) := X^fci w^X^^ starting at level j >Q is given 
by 



where the hf{Xi, . . . ,Xc) ~ J2i<ii< - <if<cl\e=i -^ie ^''"^ complete homogeneous symmetric poly- 
nomials of degree f in Xi, . . . ,Xc, denoting the small solutions of the characteristic equation (26). 

Note that the complete homogeneous symmetric polynomials satisfy the formal power series identity 

Y^h,{X,^...,X^)tf = f{,^ 

f>0 1=1 ^ 
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and thus we indeed have Tj{z) —f T{z) — 1/(1 — P{z)), for j tending to infinity. In order to fix 
the endpoint of the considered lattice paths we introduce an additional variable v encoding the steps 
S — {(1, 61), . . . , (1, 6m)}. This leads to a refined recurrence relation with respect to the refined 
characteristic polynomial 

in 

P{Xv) ^^wi,{vXf'. 

1=1 

The c small solutions of the characteristic equation l—zJ27Li (vX)'''^ = are given by the previously 
encountered (case v — I) solutions Xi(z), . . . ,Xc{z) for z in a neighborhood of zero divided by v. 
Consequently, we obtain the following result. 

Corollary 4. The generating function of meanders with steps S = {(1, 61), . . . , (1, 6^)}; weights 
n = {wi, . . . , Wjn} and characteristic polynomial P{X) := X^fci wgX^^ starting at level j > 0, where 
V marks the level of the endpoint of the path, is given by 

^ ' /=o i=l 

where the hf(Xi, . . . , Xc) are the complete homogeneous symmetric polynomials in Xi, . . . , Xc, denot- 
ing the small solutions of the characteristic equation (26). 

Remark 2. We reobtain the result of Banderier and Flajolet [1] for the enumeration of meanders by 
setting j — 0. Furthermore, the generating function of excursions, starting at level j and ending at 
level j never going below zero is obtained by extracting the coefficient [v'^]Tj{z,v). In particular, one 
can reobtain the generating function of the number of excursions starting and ending at level zero. 



8. Vicious walkers, osculating walkers, and two walks in the quarter plane 

We will show how the method of Section 3 can be used to (re-)derive some results concerning so-called 
vicious walkers and osculating walkers, and also concerning two walks confined in the quarter plane. 

8.1. Vicious walkers and osculating walkers in the lock step model. We consider three walkers 
each equipped with steps 5 = {(1,1), (1,-1)}. In the lock step model the walkers move simultaneously 
at each discrete time step such that the walkers never cross. A walker configuration is called osculating 
if the walkers never share an edge, and vicious if the walkers never meet. We also consider the model 
of up-down walkers, configurations where the first and second walker may meet and share a down 
step, and the second and third walker may also meet and share an up step. For a comprehensive 
overview of this topic we refer the reader to the work of Bousquet-Melou [5] and the references therein. 
Following [5] we consider so-called (z, j)-stars, which are non-crossing configurations of walkers with 
initial level difference 2i between the first two walkers and level difference 2j between the second and 
third walker. Let Ti j = Tij{z) = X]n>o ^"i^ j'"^" denote the generating function of the number Wn-ij 
of (i, j)-star walker configurations in the lock step model of length n. According to the step set S we 
obtain the system of recurrence relations 

T^^j = 1 + z{2Ti^j +Ti+i^j +T,^j+i + Ti^i^j+i +T,„j^i + Ti+ij^i +Ti^ij), for ij>l (27) 

with initial conditions T:!^' = T^^' = 0, i, j > for vicious walkers, tI^^ = 1 -I- z{T,^i + Ti_i,i), T^^' = 
l+z(ri j+Ti j_i), i,j > 1 for osculating walkers and conditions rj'j^' = l+z(Ti_o+7i-|_i^o+7i.i+Ti_i^i), 
Tq'^' = 1 + z{To j + Tqj+i + Ti j + Ti j-_i) for up-down walkers. 

It turns out to be beneficial to study the slightly more general system of recurrence relations 

Tij = 1 + z{wTij + Ti+ij + Tij+i + Ti-ij+i + Tij-i + + Ti^i,j), for ij > 1, (28) 

with weight w >0. For i, — > 00 we have Ti,j ^ T in the sense of formal power series, with 

T = 1 + (w + 6)zT, or equivalently T = — . 

1 — [w + ojz 
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Consequently, we set Ti j = T(l — pi j) with liiD-ij^oo Pi.j = 0. By the ansatz and the definition of 
series T we obtain the following recurrence relation for pij. 

Pl,J = Z^Wpi^j + Pi+lj + Pi,j + 1 + Pi-lj + l + Pi,j-1 + Pi+lJ-l + Pi-l,j)- 

Since we assume that linij j-^oo Pi.j — 0, we use the ansatz 

p,,, = a ■ X{zy + 13 ■ Y{zy + 7 • z{zy+\ 

with unspecified a, (3, 7 independent of i and j, assuming that \X\ < 1, |y| < 1, and \Z\ < 1 for 2; in a 
neighborhood of zero. It turns out that all series X, Y and Z satisfy the same characteristic equation, 

X = z{wX + X"^ + X + 1 + X + X^ + 1) = z{2 + {2 + w)X + 2X^), X = Y = Z. 

We obtain the following result. 

Theorem 4. A family of solutions the system of recurrence relations (28), being well defined for z 
in a neighbourhood of zero, is given in terms of the function X = X(z), defined by 

9s 1- z(2 + w)- J(l - z(2 + w)Y - 16z2 
X ^ z(2 + (2 + w)X + 2X^), with X= ^ ^ ^^-^ ^ " , 

in the following way. 

1 

1 - z(w + 6)' 

with weight w > 0, and free parameters a, (3,^ independent of i and j. 



T,,j ■ {1- a- X' ~ P ■ X^ --f- X'+^) , with T = 



Subsequently, we set w — 2 and adapt Theorem 4 to the various initial conditions. For vicious walkers 
we adapt to the initial conditions 

r£l =T(l-a-(/3 + 7)X^) = j>0, 

tI^J =T{l-l3-ia + -f)X') = i>0. 

We obtain the equations 

f3 — —7, a = —7, a ~ I, (3 — 1, 
and obtain the proper solution a = (3 = —7 = 1. For osculating walkers we adapt to the initial 
conditions T^^' = 1 + z{Ti^^ + Tij^i) and tJq' = 1 + z{T,^i + Ti„i_i), leading to the system of 
equations 

T(l - a - (/3 + -/)X^)) = 1 + zT(l -aX - I3X^ - lX^+^) + zT{l - aX - (3X^-^ - -^X^) j > 0, 

T(l - /3 - (a + 7)X') = 1 + zT(l - aX' ~ (3X ~ -fX'+^) + zT{l - aX*"^ ~ (3X - 7X') i > 0. 
We obtain the equations 

r(l - a) = 1 + 2zT(l - aX), T(l - /3) = 1 + 2zT(l - (3X), 

(a + 7) = z{a{l + ^)+ 7(1 + X)),{I3 + 7) = z(/3(l + 1) + 7(1 + X)), 

and easily obtain the proper solution a — /3 — 3X/{1 + 2X), 7 ~ SX/{2 + X), using the relations 
z = X/{2{1 + X)'^) and T — 1/(1 — 8z). For up-down-walkers we obtain the equations 

T(l - a) = 1 + 2zT{2 -a- aX), T{1 - (3) = 1 + 2zT{2 - (3 ~ [3X), 
(a 7) = z{a{2 + X+^) + 2^{l + X)), (/3 7) = z{l3{2 + X + ^) + 2^{l + X)), 
and easily obtain the solutions a = (3 = 2X/ (1 -|- X) and 7 = —X. 

Corollary 5 (Bousquet-Melou [5], Gessel). The length generating function of three vicious walkers, 
three osculating walkers and three up-down walkers in (i,j)-star configuration is given in terms of 
series X , given by 

1 - 4z - VI -8z 



X — 2z{l + X) , or equivalently X 



4:Z 
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in the following way. For vicious walkers 



rr ( ^ rr.(, / ^ ^ 2{l + X) - u{l + wX) + A 
r,, (z, u,w)^T{^l- a{X + X=) +^ ^,^.x)- uXil + w) ' 



for osculating walkers 

■''J ^ 1 + 2X 1 + 2X 2 + X ' l-Sz^ 1 + 2X 1 + 2X 2 + X 

and for up-down walkers 

''^ ^ l + X l + X ' l-8z^ l + X l + X ' 

One can easily obtain a refinement of the results above by counting the number of osculations with 
variable u and the number of shared up and down edges by variable w, which leads to a refined of 
Theorem 4. The results presented above correspond to the special cases (u, w) = (0, 0) , {u,w) = (1, 0) 
and = (1, 1), respectively. 

Corollary 6. The length generating function generating function of stars where u counts the 

number of osculations and w the number of shared up-down edges after osculations, is given by 

2(1 + X) - uil + wX) 
'2{l + X)-uX{l + w) 

with a given by the following expression. 

(1 + X)2 -u(l-X + X^ + wX) 
(l + X)2-uX(w + X) 

8.2. Vicious walkers and osculating walkers in the random turn model. The random turn 
model of walkers is similar to the lock step model but at each discrete time step only one of the 
walkers is allowed to move. We consider three vicious walkers and osculating walkers with Dyck steps 
S = {(1, 1), (1, -1)} or Motzkin steps S = {(1, 1), (1, -1), (1,0)}. Let T,,, = T,,,{z) = J2n>o Wn■,^,,z^ 
denote the generating function of the total number of Wn-ij the number of walker configurations of 
length n, with initial level difference i between the first two walkers and level difference j between 
the second and third walker. For the step set iS = {(1, 1), (1, —1)} we obtain the system of recurrence 
relations 

Tij = 1 + z(Ti+i,j + Ti,j+i + Ti_ij+i + Tij^i + Tj+ij-i + Ti-ij), for i,j >1 (29) 

with initial conditions t[^^ = T^^j = 0, i,j > for vicious walkers, and Ti._i = T^i.j = for 
osculating walkers. For the step set S — {(1, 1), (1, — 1), (1, 0)} we obtain the system of recurrence 
relations 

ri,j = 1 + z(3rj+ij + T:j.j+i + T:,_ij+i + rj,j_i + Ti+ij_i + Tj-ij), for i,j>l (30) 

with initial conditions tIq' = T^^j = 0, «,j > for vicious walkers, and = tI'^'^ = for 

osculating walkers. Using Theorem 4 we obtain the following results. 

Corollary 7. The length generating function of three vicious walkers and osculating walkers in the 
random turn model with Dyck steps S — {(1, 1), (1, — 1)} is given in terms of series X, defined by 
X — 2z{\ + X + X'^). For vicious walkers we obtain tJ^' = ' (1 ~ ^')(1 ^ X^). The generating 
function of osculating walkers is given by Tf^^{z) = T!j'^[ In particular, we obtain 

Corollary 8. The length generating function of three vicious walkers and osculating walkers in the 
random turn model with Motzkin steps S = {(1, 1), (1, —1), (1, 0)} is given in terms of series X, defined 
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by X = z{2 + 5X + 2X'^). For vicious walkers we obtain i;'^' = • (1 - X'){1 - X^), and for 

osculating walkers we get tJ'j'' — jzx^ ' (1 ^ X^^^)(l — X^^^). In particular, we obtain 



8.3. Two families of walks confined in the quarter plane. Recently, Bousquet-Melou and 
Mishna [7] presented a systematic approach to the problem of counting walks starting at the origin 
(0, 0) being confined to the quarter plane NoxNq. In particular, for step sets Si = {(—1,0), (0, 1), (1, —1)} 
and ^2 = {(—1,0), (0, 1), (1, 0), (0,-1), (—1, 1), (1,-1)} they derived, amongst many other results, the 
generating function of walks of length n starting at the origin (0, 0), the result involves the generating 
function of the Motzkin numbers. For the two models Si and ^2 we will derive the generating function 
of the number of paths starting at confined to the quarter plane No x Nq using the approach of 
Bouttier, Di Francesco and Guitter; note that the subsequently derived results can also be obtained, 
in a more systematic manner, using the kernel method [7]. Let Ti,j = Tij{z) denote the generating 
function of the total number of paths starting at point (i, j) G No x No confined to the quarter plane. 
Concerning the first step set Si — {(—1, 0), (0, 1), (1, —1)} we obtain the system of recurrence relations 
i0TT,,,^Tg'\z), 

Tij(z) = 1 + z(T'i_ij + Tij+i + Ti+ij.i), fori,j>0, (31) 

with initial values T_i j — Ti _i = for i,j > 0. We set T^ j — T(l — pij) with linii j^^pi j — 0, 
and series T given as the solution of the limiting equation T = 1 + 3zT. By the ansatz we obtain the 
following recurrence relation for pij. 

Pi,j = ^(Pi-ij + Pi,j+i + Pi+i,j)- 

Since we assume that limi j^oo Pi,j — 0, we use again the ansatz 

p,,, = a ■ X{zy + 13 ■ Y{zy + 7 • z{zy+\ 

assuming that \X\ < 1, \Y\ < 1, and |Z| < 1 for 2; in a neighborhood of zero. It turns out that all 
three series satisfy the same characteristic equation, 

X^z{l + X + X^), or equivalently X = ^ ' ^ ' + ^^)(l-3^) ^ ^.^^ X = Y = Z. 

Note that X is the generating function of the Motzkin numbers. Hence, we obtain the family of 
solutions 

1 

1 ~3z' 

Adapting to the initial conditions T^ij — = 0, for i,j > 0, gives a = (3 = X and 7 = —X^. 

For the second step set ^2 = {(— 1, 0), (0, 1), (1, 0), (0, — 1), (— 1, 1), (1, — 1)} we obtain the system of 

rc 1 

recurrence r elations for Tij = Tlj'{z), 

Ti^-j 1 + 2:(Ti_ij + Ti+i,j + Tij+i + Tij^i + Tj-ij+i + Ti+ij-i), for i,j > 0, 

with initial values T^i.j = Ti__i = for i,j > 0. This model is equivalent to the random turn model 
of three osculating Dyck walkers, and already solved in Corollary 7. 

r o ] 

Theorem 5. The length generating functions T-j{z) of paths confined to the quarter plane starting 
at {i,j) e No X No with step set Si = {(—1, 0), (0, 1), (1,-1)} is given by 

Tlfiz) = ^^(1 - Xl+') (1 - Xi+'), with Xl = z{l + Xi+ Xf). 

The length generating functions rj^^'(z) of paths confined to the quarter plane starting at {i,i) G 
No X No with step set S2 = {(^1, 0), (0, 1), (1, 0), (0, —1), (—1, 1), (1, ~1)} equals the length generating 
function of three osculating Dyck walkers in the random turn model, 

7;5^i (z) = rg^i (2z) = (1 - x^+i) (1 - xi+') , with X2 = 2z(i + X2 + xl). 



= ^(1 - a • - /3 • X^' - 7 • X'+^ 
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Conclusion 

We have shown that the "asymptotic series method" of Bouttier, Di Francesco and Quitter can be used 
to study several families of embedded trees. Moreover, we used the method to study simple families 
of lattice path which can be considered as a degenerated family of embedded trees. Furthermore, we 
presented some other problems which can be solved using this method together with a suitable ansatz. 
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